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ABSTRACT 



A thorough search for large scale anisotropies in the distribution of arrival directions 
of cosmic rays detected above 10 18 eV at the Pierre Auger Observatory is reported. For 
the first time, these large scale anisotropy searches are performed as a function of both 
the right ascension and the declination and expressed in terms of dipole and quadrupole 
moments. Within the systematic uncertainties, no significant deviation from isotropy 
is revealed. Upper limits on dipole and quadrupole amplitudes are derived under the 
hypothesis that any cosmic ray anisotropy is dominated by such moments in this energy 
range. These upper limits provide constraints on the production of cosmic rays above 
10 18 eV, since they allow us to challenge an origin from stationary galactic sources 
densely distributed in the galactic disk and emitting predominantly light particles in all 
directions. 



Subject headings: astroparticle physics — cosmic rays 



The large scale distribution of arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) 
as a function of the energy is a key observable to provide further understanding of their origin. 

Above ~ 0.25 EeV, the most stringent bounds ever obtained on the dipole com ponent in the equato- 

rial plane were recently reported, being below 2% at 99% C.L. for EeV energies (jAuger Collaboration 201 lal ) 
Such a sensitivity provides some constraints upon scenarios in which dipolar anisotropies could be 
imprinted in the distribution of arrival directions as the result of the es c ape of UHECRs from the 
Galaxy up to the ankle energy (jPtuskin et al. 19931 : ICandia et al. 20031 : iGiacinti et al. 20121). On 
the other hand, if UHECRs above 1 EeV have already a predomina nt extragalactic origin (IHillas 19671 : 
Blumenthal 19701 : iBerezinskv et al. 20061 : iBerezinskv et al. 2004 ) , their angular distribution is ex- 
pected to be isotropic to a high level. Thus, the study of large scale anisotropies at EeV energies 
would help in establishing whether the origin of UHECRs is galactic or extragalactic in this energy 
range. 

The upper limits aforementioned are based on first harmonic analyses of the right ascension 
distributions in several energy ranges. The analyses benefit from the almost uniform directional 
exposure in right ascension of any ground based observatory operating with high duty cycle, but 
are not sensitive to a dipole component along the Earth rotation axis. In contrast, using the large 
amount of data collected by the surface detector array of the Pierre Auger Observatory, we report 
in this letter on searches for dipole and quadrupole patterns significantly standing out above the 
background noise whose components are functions of both the right ascension and the de clination 
(a detailed description of the present analysis can be found in (jAuger Collaboration 20121 ) ). 
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The Pierre Auger Observatory is located in Malargiie, Argentina, at mean latitude 35.2° S, 
mean longitude 69.5° W and mean altitude 1400 meters above sea level. It exploits two available 
techniques to detect extensive air showers initiated by UHECRs : a Surface Detector (SD) array 
and a Fluorescence Detector (FD). The SD array consists of 1660 water-Cherenkov detectors laid 
out over about 3000 km 2 on a triangular grid with 1.5 km spacing, sensitive to the light emitted 
in their volume by the secondary particles of the showers. At the perimeter of this array, the 
atmosphere is overlooked on dark nights by 27 optical telescopes grouped in 5 buildings. These 
telescopes record the number of secondary charged particles in the air shower as a function of depth 
in the atmosphere by measuring the amount of nitrogen fluorescence caused by those particles along 
the track of the shower. At the lowest energies observed, the angular resolution of the SD is about 
2.2°, and reaches ~ 1° at the highest energies. This is sufficient to perform searches for large- 
scale anisotropics. The statistical fluctuation in energy measurement amounts to about 15%, while 
the a bsolute energy scale is giy en by the FD measurements and has a systematic uncertainty of 



22% (jAuger Collaboration 2008 ). 



In the analyses presented in this letter, the data set consists of events recorded by the SD 
array from 1 January 2004 to 31 December 2011, with zenith angles less than 55°. To ensure good 
reconstruction, an event is accepted only if all six nearest neighbours of the water-Cherenkov detec - 



tor with the highest signal were operational at the time of the event (jAuger Collaboration 2010al ) . 
Based on this fiducial cut, any active water-Cherenkov detector with six active neighbours de- 
fines an active elemental cell. In these conditions, and a bove the energy at which the detection 



efficiency saturates, 3 EeV (jAuger Collaboration 2010a! ) . the total exposure of the SD array is 
23,520 km 2 yr sr. 

Due to the steepness of the energy spectrum, any mild bias in the estimate of the shower energy 
with time or zenith angle can lead to significant distortions of the event counting rate above a given 
energy. It is thus critical to control the energy estimate in searching for anisotropies. The procedure 
followed to obtain an unbiased estimate of the sho wer energy consists in correc ting measurements of 



showe r signals for the influences of weather effects (jAuger Collaboration 20091 ) an d the geomag netic 



the 



field (jAuger Collaboration 201 lbl ). Using the constant intensity cut method (jrlersil 1961 
shower signal is then converted to the value that would have been expected had the shower arrived 
at a zenith angle 38° . This reference shower signal is finally converted into energy using a calibration 
curve based on hybrid events measured simultaneously by the S D array and FD telescopes , since 



the latter can provide a calorimetric measurement of the energy ([Auger Collaboration 2008! ) 



In searching for anisotropies, it is also critical to know accurately the effective time-integrated 
collecting area for a flux from each direction of the sky, or in other words, the directional expo- 
sure oj of the Observatory. For each elemental cell, this is obtained through the integration over 
Local Sidereal Time (LST) a of iW(a°) x a ce ii(0) x e(9,(p,E), with x^\a°) the total operational 
time of the cell (i) at LST a , a ce ii(#) = 1.95 cosfl km 2 the geom etric aperture of each elemen- 



tal cell under incidence zenith angle 9 (jAuger Collaboration 2010al ). and e(6,ip,E) the detection 



efficiency under incidence zenith angle 9 and azimuth angle cp at energy E. In the same way as 
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in ( Auger Collaboration 201 lal ). the small modulation of the exposure in a due to the variations 
of scW can be accounted for by re- weighting the events with the number of elemental cells at the 
LST of each event k, AN ce \\(a®). Since both 8 and ip depend only on the difference a — a , the 
integration over a can then be substituted for an integration over the hour angle a' = a — a so 
that the directional exposure actually does not depend on right ascension when the zW are assumed 
to be independent of the LST : 

"coll ,-2Ah 

lj(6, E)=J2 x« / da' a ce ii(0(a', 5)) e(0(a', 5), tp(a/,S),E). (1) 
i=i Jo 

The zenithal depende nce of the detection efficien cy e(6, tp, E) can be obtained directly from the data 



in an empirical way (jAuger Collaboration 20121 ) . Additional effects have an impact on u), such as 
the azimuthal dependence of the efficiency due to geomagnetic effects, the corrections to both the 
geometric aperture of each elemental cell and the detection efficiency due to the tilt of the array, 
and the corrections due to the spatial extension of the arr ay. Accounting for all thes e effects, the 



resulting dependence of uj on declination can be found in (jAuger Collaboration 20121 ) . For a wide 
range of declinations between ~ —89° and ~ —20°, the directional exposure is ~ 2,500 km 2 yr at 
1 EeV, and ~ 3, 500 km 2 yr for any energy above full efficiency. Then, at higher declinations, it 
smoothly falls to zero, with no exposure above 20° declination. 

The detection of significant dipole or quadrupole moments above EeV energies would be of 
considerable interest. Dipole and quadrupole patterns are encoded in the low order a% m and a2m 
coefficients of the multipolar expansion of any angular distribution over the sphere <£>(n) : 

$(n) = ^ ^m^m(n), (2) 

where n denotes a unit vector taken in equatorial coordinates. Due to the non-uniform and in- 
complete coverage of the sky at the Pierre Auger Observatory, the estimated coefficients a,£ m are 
determined in a two-step procedure. First, from any event set with arrival directions {ni, ...,iin} 
recorded at LST {a°, a^}, the multipolar coefficients of the angular distribution coupled to the 
exposure function are estimated through : 

N v i \ 

AiV ce n(a!j!) corrects for the slightly non- uniform directional exposure in right ascension. Then, 
assuming that the multipolar expansion of the angular distribution <£(n) is bounded to £ max , the 
first bi m coefficients with £ < £ max are related to the non-vanishing a^ m through : 



W&fvem?, (4) 



p=0 m '= 
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Fig. 1. — Reconstructed amplitude of the dipole as a function of the energy. The dotted line stands for the 
99% C.L. upper bounds on the amplitudes that would result from fluctuations of an isotropic distribution. 



where the matrix K is entirely determined by the directional exposure : 

[K]f™' = [ dO w(n) Y im (n) *W(n). (5) 
J An 

Inverting Eqn. 0] allows us to recover t he underlying a^ m , with a resolution proportional 



to ([K^ 1 ]^™ aoo) ' 5 ( Billoir Deligny 20081 ). As a consequence of the incomplete coverage of the 
sky, this resolution deteriorates by a factor larger than 2 each time £ max is incremented by 1. With 
our present statistics, this prevents the recovery of each coefficient with good accuracy as soon as 
f max > 3, which is why we restrict ourselves to dipole and quadrupole searches. 

We first assume that the angular distribution of cosmic rays is modulated by a pure dipole 
and parameterise the intensity <I>(n) in any direction as : 

$(n) = ^|l|rdTiJ, (6) 

where d denotes the dipole unit vector. The dipole pattern is here fully characterised by a dec- 
lination 5d, a right ascension a^, and an amplitude r corresponding to the maximal anisotropy 
contrast : r = (<& ma x — $min)/($max + ^mm)- The estimation of these three coefficients is straight- 
forward from the estimated spherical harmonic coefficients a\ m - The reconstructed amplitudes r 
are shown in Fig. [1] as a function of the energy. The 99% C.L. upper bounds on the amplitudes 
that would result from fluctuations of an isotropic distribution are indicated by the dotted line. 
One can see that within the statistical uncertainties, there is no evidence of any significant signal. 
In Fig. [21 the corresponding directions are shown in orthographic projection with the associated 
uncertainties, as a function of the energy. Both angles are expected to be randomly distributed 



-10- 



*x=0° 



o 





8=30° 




a=30° 
8=0° 








8=-30° 








□ 

8=-60° . 











a=60° 



O 1<E[EeV]<2 
□ 2<E[EeV]<4 
O 4<E[EeV]<8 
E[EeV]>8 



Fig. 2. — Reconstructed declination and right-ascension of the dipole with corresponding uncertainties, as 
a function of the energy, in orthographic projection. 



in the case of independent samples whose parent distribution is isotropic. It is thus interesting to 
note that all reconstructed declinations are in the equatorial southern hemisphere, and to note also 
the intriguing smooth alignment of the phases in right ascensio n as a function of the energy . In our 
previous report on first harmonic analysis in right ascension ([Auger Collaboration 201 lal ). we al- 
ready pointed out this alignment, and stressed that such a consistency of phases in adjacent energy 
intervals is expected with smaller number of events than the detection of amplitudes standing-out 
significantly above the background noise in the case of a real underlying anisotropy. This motivated 
us to design a prescription aimed at establishing at 99% C.L. whether this consistency in phases 
is real, using the exact same analysis as the one reported in (j Auger Collaboration 201 lal ). The 
prescribed test will end o nce the total exposure since 25 June 2011 reaches 21,000 km 2 yr sr. The 
smooth fit to the data of (lAuger Collaboration 201 lal ) is shown as a dashed line in FigO restricted 
to the energy range considered here. Though the phase between 4 and 8 EeV is poorly determined 
due to the corresponding direction in declination pointing close to the equatorial south pole, it is 
noteworthy that a consistent smooth behaviour is observed using the analysis presented here and 
applied to a data set containing two additional years of data. 

Assuming now that the angular distribution of cosmic rays is modulated by a dipole and a 
quadrupole, the intensity <I>(n) can be parameterised in any direction n as : 



$(n) 



47T 



1 + r d • n + A + (q H 



n 



+ A (qo -n) 2 + A_(q_ • n) 
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Fig. 3. — Reconstructed right asce nsion of the dipole as a function of the energy. The smooth fit to the 
data of ([Auger Collaboration 201 lal) is shown as the dashed line (see text). 



with the constraint A+ + A- + Ao = 0. It is convenient to define the quadrupole amplitude f3 = 
(A+ — A-)/ (2 + A+ + A_), which provides a measure of the maximal quadrupolar contrast in the 
absence of a dipole. Hence, any quadrupolar pattern can be fully described by two amplitudes 
(/?, A-f) and three angles : (<5+,a+) which define the orientation of q+ and (a_) which defines the 
direction of q in the orthogonal plane to q+. The third eigenvector qo is orthogonal to q+ and 
q . The estimated amplitudes A + and j3 are shown in Fig. [J] as functions of the energy. In the same 
way as for dipole amplitudes, the 99% C.L. upper bounds on the quadrupole amplitude that could 
result from fluctuations of an isotropic distribution are indicated by the dashed lines. Throughout 
the energy range, there is no evidence for anisotropy. 

There are small uncertainties in correcting the estimator of the energy for weather and geomag- 
netic effects, and these propagate into systematic uncertainties in the measured anisotropy parame- 
ters. As well, anisotropy parameters may be altered in a systematic way by energy dependence of th e 



attenuation curve. All these systematic effects have been quantified ([Auger Collaboration 20121 ) 
They do not change significantly the results presented here. 



From these analyses, upper limits on dipole and quadrupole amplitudes can be derived at yyyo 
C.L.. They are shown in Fig. [5] for the dipole amplitudes, accounting for the systematic uncer- 
tainties. We illustrate now their astrophysical interest by calculating the amplitudes of anisotropy 
expected in a toy scenario in which sources of EeV-cosmic rays are stationary, densely and uniformly 
distributed in the galactic disk, and emit particles in all directions. 

Both the strength and the structure of the magnetic field in the Galaxy, known only approx- 
imately, play a crucial role in the propagation of cosmic rays. The field is thought to contain a 
large scale regular component and a small scale turbulent one, both having a local strength of a 
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Fig. 4. — Amplitudes of the quadrupolar moment as a function of the energy using a multipolar reconstruc- 
tion up to £ max — 2. The dotted lines stand for the 99% C.L. upper bounds on the amplitudes that could 
result from fluctuations of an isotropic distribution. 



few microgauss (see e.g. (jBeck 200 ll )). While the turbulent component dominates in strength by 
a factor of a few, the regular component imprints dominant drift motions as soon as the Larmor 
radius of cosmic rays is larger than the maximal scale of the turbulences (thought to be in the 
range 10-100 pc). We adopt here a recent parameterisation of the regular component obtained 
by fitting model field geometries to Faraday rotation measures of extragalactic radio sources and 
polaris ed synchrotron emissi on (BSS-model, with anti-symmetric halo with respect to the galactic 
plane) (jPshirkov et al. 201ll ). In addition to the regular component, a turbulent field is generated 
according to a Kolmogorov power spectrum and is pre-computed on a three dimensional grid peri- 
odically repeated in space. The size of the grid is selected to match the maximal scale of turbulences 
(taken here as 100 pc), and the strength of the turbulent component is taken as three times the 
strength of the regular one. To describe the propagation of cosmic rays with energies E > 1 EeV in 
such a magnetic field, the direct integration of trajectories is the most appropriate tool. To obtain 
the anisotropy of cosmic rays emitted from sources uniformly distributed in a cylinder with a radius 
of 20 kpc from the galactic cent re and with a height of ± 100 pc, we adopt a method first proposed 
in (jThielheim &: Langhoff 19681 ). It consists in back-tracking anti-particles with random directions 
from the Earth to outside the Galaxy. Each test particle probes the total luminosity along the path 
of propagation from each direction as seen from the Earth. For stationary sources emitting cosmic 
rays in all directions, the expected flux in the initial sampled direction is proportional to the time 
spent by each test particle in the source region. 

The amplitudes of anisotropy obviously depend on the rigidity E/Z of the cosmic rays, with Z 
the electric charge of the particles. Since we only aim at illustrating the upper limits, we consider 
two extreme single primaries : protons and iron nuclei. The calculation of anisotropy amplitudes 
for single primaries is useful to probe the allowed contribution of each primary as a function of the 
energy. 
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Fig. 5. — 99% C.L. upper limits on dipole and quadrupole amplitudes as a function of the energy. Some 
generic anisotropy expectations from stationary galactic sources distributed in the disk are also shown, for 
various assumptions on the cosmic ray composition. The fluctuations of the amplitudes due to the stochastic 
nature of the turbulent component of the magnetic field are sampled from different simulation data sets and 
are shown by the bands. 



The dipole and quadrupole amplitudes obtained for several energy values covering the range 
1 < E/~EeV < 20 are shown in Fig. [5j To probe unambiguously amplitudes down to the percent 
level, it is necessary to generate simulated event sets with at least ~ 5 10 5 test particles. Such a 
number of simulated events allows us to shrink statistical uncertainties on amplitudes at the 0.5% 
level. Meanwhile, there is an intrinsic variance in the model for each anisotropy parameter due to 
the stochastic nature of the turbulent component of the magnetic field. This variance is estimated 
through the simulation of 20 sets of 5 10 5 test particles, where the configuration of the turbulent 
component is frozen in each set. The RMS of the amplitudes sampled in this way is shown by the 
bands in Fig. 

The resulting amplitudes for protons largely stand above the allowed limits. Consequently, 
unless the strength of the magnetic field is much higher than in the picture used here, the up- 
per limits derived in this analysis exclude that the light component of cosmic rays comes from 
galactic stationary sources densely distributed in the galactic disk and emitting in all directions. 
To respect the dipole limits below the ankle energy, the fraction of protons should not exceed 
~ 10% of the cosmic ray composition. This is particularly interesting in the view of the indi- 
cations for the presence of a light component around 1 EeV from shower depth maximum mea- 



surements (lAuger Collaboration 2010bl : lAbbasi et al. 20101 : lJui et al. 20 111 ) , though firm interpre- 



tations of these measurements in terms of the atomic mass still suffer from some ambiguity due to 
the uncertain hadronic interaction models used to describe the shower developments. On the other 
hand, if the cosmic ray composition around 1 EeV results from a mixture containing heavy elements 
of galactic origin and light elements of extragalactic origin, upper limits can be respected. This 
is because large scale anisotropy amplitudes below the percent level are expected for extragalactic 
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cosmic rays, du e to the motion of the Galaxy relative to a poss ibly stationary extragalactic cosmic 
ray rest frame (jKachelriess Serpico 2000 ; lHarari et al. 20101 ) . 



Future measurements of composition below 1 EeV will come from the l ow energy extension 

Combining these 



HEAT now available at the Pierre Auger Observatory (IMathes et al. 2011 



measurements with large scale anisotropy ones will then allow us to further understand the origin 
of cosmic rays at energies less than 4 EeV. 
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